Topological Order with a Twist: Ising Anyons from an Abelian Model 
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Anyon models can be symmetric under some permutations of their topological charges. One 
can then conceive topological defects that, under monodromy, transform anyons according to a 
symmetry. We study the realization of such defects in the toric code model, showing that a process 
where defects are braided and fused has the same outcome as if they were Ising anyons. These ideas 
can also be applied in the context of topological codes. 
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Particle statistics are particularly rich in two spatial di- 
mensions, where beyond the usual fermions and bosons 
there exist more generally anyons (see [l[ for a compi- 
lation of the basic references). Anyonic statistics are 
complex enough to give rise to the notion of topologi- 
cal quantum computation (TQC) 0-01: where compu- 
tations are carried out by braiding and fusing anyons, 
see Fig. [Tfa-c). The nonlocal encoding of quantum in- 
formation on fusion channels and the topological nature 
of braiding makes TQC naturally robust against local 
perturbations, providing a complement to fault-tolerant 
quantum computation [^, 0] . 

In condensed matter, anyons emerge as excitations in 
systems that exhibit topological order 0. A possible 
way to obtain these exotic phases is by engineering suit- 
able Hamiltonians on lattice spin systems iiSli. In- 
deed, implementations on optical lattices have been pro- 



posed [12|. Unfortunately, the anyon models that ap- 



pear in simple models are not computationally powerful. 
In this paper we address an strategy to recover compu- 
tationally interesting anyon-like behavior from systems 
with very simple anyonic statistics. 

Our starting point are the symmetries that anyons may 
exhibit. Anyon models have three main ingredients: (i) 
a set of labels that identify the superselection sectors or 
topological charges, (ii) fusion/splitting rules that dictate 
the charges of composite systems, and (iii) braiding rules 
that dictate the effect of particle exchanges. A symme- 
try is a label permutation that leaves braiding and fusion 
rules unchanged — for a recent survey, see [13| — . Given a 
symmetry s, we can imagine cutting the system along an 



'»»»»/ 



(a) 



(b) 



(d) 



FIG. 1: Anyon processes: time flows upwards. (a,b) Two 
topologically distinct ways to exchange anyons. (c) A fusion 
of two anyons. (d) Two twists (crosses) connected by a line 
(dotted) across which e charges become m charges. 



open curve, as in Fig.[ljd), and then gluing it again "up 
to s" . Ideally the location of the cut itself is unphysical, 
only its endpoints have a measurable effect. In partic- 
ular, transporting an anyon around one end of the line 
changes the charge of the anyon according to the action 
of s. Our aim is to explore to which extent these topolog- 
ical defects, that we call twists for short, can be "treated 
as anyons" and used in TQC. Twists are being indepen- 
dently studied by Kong and Kitaev 1J|. An interesting 
precedent are the Alice strings appearing in some gauge 
models [TiJ, which can cause charge conjugation under 
monodromy, whereas the twists that we will discuss here 
exchange electric and magnetic charges. 

Rather than trying a general, abstract approach, we 
will focus on a well-known spin model, the toric code 
model, and address twists constructively. In this model 
anyons have no computational power, but we will show 
that twists behave as Ising anyons 16], which are compu- 
tationally interesting. In fact, they do not directly allow 
universal computation, but there exist strategies to over- 
come this difficulty [l7l - [l9j |. In [2^, Wootton et al. also 
try to mimic the non-abelian behavior in an abelian sys- 
tem, using an entirely different approach and philosophy. 

We remark that, although the discussion will mainly 
be in terms of topological order, it has direct application 
in the closely related context of topological codes (2l| - [2^ . 

Anyon models — Anyon models are mathematically 
characterized by modular tensor categories, but we will 
not need such generalities (for an introduction, see for 
example [1^). Instead, we will illustrate the content of 
anyon models with an example: Ising anyons. 

The first element of an anyon model is a set of la- 
bels that identify the superselection sectors or topological 
charges of the model. For Ising anyons there are three: 
1, a and "tp. Any given anyon carries such a charge, which 
cannot be changed locally. We can also attach a charge 
to a set of anyons or a to a given region. A region without 
anyons has trivial charge 1. 

Next we need a set of fusion rules that specify the pos- 
sible values of the total charge in a composite system. In 
terms of anyon processes, fusion rules specify the possible 
outcomes of the fusion of two anyons, see Fig. [TJc). For 



Ising anyons fusion rules take the form 

a X a = 1 + ip, a X tp = (7, ip x tp 



(1) 



That is, a pair of ct-s may fuse into the vacuum or produce 
a ■(/', a (T and a ip always fuse into a and two ip-s into the 
vacuum. Fusion rules are commutative and 1 x a = a. 

When two a anyons are far apart, their total charge, 
which might be 1 or -0, becomes a non-local degree of 
freedom. This is indeed an example of a topologically 
protected qubit, since there are two possible global states. 
We can measure this qubit in the charge basis by fusing 
the two cr-s and checking the output. In general for any 
set of anyons with given charges there is a fusion space 
that describes the non-local degrees of freedom related to 
fusion outcomes. For example, for 2n as with indefinite 
total charge the fusion space has dimension 2". 

Braiding operations as those in Fig. [lja,b) act on the 
fusion space in a topologically protected way. This action 
is in general described by braiding rules, but in the case 
of Ising anyons it is possible to characterize braiding with 
Majorana operators. In particular, for 2n as we need 2n 
majorana operators These are self-adjoint operators 
that satisfy cjCk+CkCj — 25 jk and act on , as needed. 
The total charge of the j-th and j + 1-th anyon is given 
by the operator — zcjCj+i, the eigenvalues -fl and —1 
corresponding to the total charge 1 and ip, respectively. 
Under the braiding of Fig. [TJa) the operators evolve as 
follows: 
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FIG. 2: A square lattice with spins living at vertices, (a) Pla- 
quette operators are products of four Pauli operators. String 
operator are products of Pauli operators on their edges, (b) 
A dislocation in the geometry of the Hamiltonian produced 
by shifting plaquettes. In the pentagon one can introduce the 
indicated plaquette operator, which commutes with the rest. 



A remarkable property of the rules ([31 S]) is their invari- 
ance under the exchange of e and to. Indeed, there exists 
a self-equivalence of the corresponding modular tensor 
category that exchanges e and m ^1^. It is then con- 
ceivable a construction as that in Fig. [IJd) , where there 
is a line in the system across which e and m charges are 
exchanged. The rest of this paper is devoted to explicitly 
realize this in a spin model and explore its consequences. 

Toric code model — This is a spin-1/2 Hamiltonian 
model, with spins forming a square lattice. The inter- 
actions are four-body, with each plaquette in the lattice 
contributing a Pauli product term to the Hamiltonian: 
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(2) 



Ak := ala 



(5) 



where j ^ k =^ j + 1. This describes braiding up to 
process dependent phases that we will not need. 

The quantum gates obtained from the braiding of a 
anyons are not universal. Indeed, they are Clifford gates. 
Yet, they can be complemented with physically plausi- 
ble, topologically unprotected, noisy operations to get 
universal quantum computation 17]. 

Symmetries — We will now illustrate the notion of a 
symmetry in an anyon model with the model technically 
known as the quantum double of the group Z2. It has 
four charges: 1, e, m and e. The fusion rules are 

exm = e, exe = m, m x e — e, 

exe = mxm — exe — 1. (3) 

These are abelian fusion rules: the result of a fusion has 
always a definite outcome. Thus, the fusion space is triv- 
ial and the effect of braiding a charge a with a charge b 
as in Fig.fTfa) can be captured in a phase Rab- The Rab 
phases have in general no measurable effect, but we can 
build the following invariants, in principle measurable: 

Ree — Rmm — I5 -^ee li 

Rem-^me — Ree-^ee — Rnie-^em 1- (4) 



Here k = {i,j) indexes the spins and i := (1,0), j := 
(0, 1). We have chosen a uniform coupling J — 1; taking 
instead a different non-zero coupling for each term would 
not change the physics that interest us. The Hamilto- 
nian ([S|) might be hard to engineer, but it can be ob- 
tained effectively from a two-body spin-1/2 model (8|] for 
which there exist experimental proposals 12|. The toric 
code model was originally introduced by Kitaev 0], but 
the more symmetric form ^ is due to Wen This 
symmetry was further studied in [26 1 . 

The ground subspace of ([5]) is described by the con- 
ditions Ak = 1. Excitations are localized and gapped: 
a plaquette k is excited if A^ = —1, in which case we 
say that it holds a quasiparticle. These quasiparticles are 
anyons, described by the quantum double of Z2 (|3l4p . To 
label quasiparticles with their topological charge, we first 
have to label plaquettes with two 'colors' as in a chess- 
board lattice, see Fig.[2Ia). Then we can attach a charge 
e (charge to) to quasiparticles living at dark (light) pla- 
quettes. Notice that the exchange of e and m labels is 
trivially a symmetry at the Hamiltonian level, since the 
choice of dark/light plaquettes is entirely arbitrary. 
String operators — Consider the action of a^. It fiips 
Thus e and to charges are bosons and e charges fermions. the state of the plaquettes k — i and fc— j, either hopping a 



quasiparticle or creating/ annihilating a pair of quasipar- 
ticles. We can thus visuahze as a segment t connecting 
the centers of these two plaquettes and write St cr^- 
Similarly, we relate to a segment t' connecting the centers 
of /c and fc — i— j the operator Sf := cr^- Next we consider 
a string 7, an ordered collection of segments and 
build the string operator S-y St^ ■ ■ • S'tj^ . Given an open 
string 7, with two different plaquettes k and k' as end- 
points, flips the eigenvalues of Ak and Ak'. Thus, S-y 
represents a process where an anyon is moved along the 
string, plus suitable creation/ annihilation events on the 
endpoints. We can distinguish between dark and light 
strings, depending on whether they connect light or dark 
plaquettes; dark (light) strings transport e anyons (m 
anyons). String operators have two important proper- 
ties: (i) a dark string 7-° and a light string 7^ anticom- 
mute when they cross and odd number of times, and (ii) 
if Vr is the subspace of states with trivial charge in a 
given region R and two strings 7 and 7' enclose R, then 
SjIvr = Sj'lvit- That is, strings operators can be freely 
deformed as long as they do not go over an excitation. 

The algebra of string operators alone is enough to com- 
pute the anyonic nature of the model [1, . Here we are 
more interested on the fact that among string operators 
there are constants of motion related to the topological 
charge of given regions. Notice that closed string oper- 
ators commute with all plaquette operators; they rep- 
resent vacuum to vacuum processes. Consider a closed 
dark string 7-° and a closed light string 7^ enclosing a 
given region. Then S^d = 1 if the total charge is 1 or 
e, S.yD = — 1 otherwise. Similarly, S^l = 1 if the to- 
tal charge in the region is 1 or m, S^n = — 1 otherwise. 
Together, S^d and S^l distinguish the four charges. 

Twists — To construct a line L across which e and 
m are exchanged, as in Fig. [Ijc), we have to introduce 
dislocations in the lattice: along L, plaquettes are shifted 
so that the coloring does not match, as in Fig.[2]Jb). The 
geometrical change implies a change in the Hamiltonian, 
and we indicate in the figure the new plaquette operators. 
At the ends of the line we find pentagonal plaquettes. 
We can either add a pentagonal plaquette operator, as 
indicated in the figure, or allow a localized gapless mode. 
Turning points of L can be dealt with similarly. 

The topologically interesting behavior occurs at the 
ends of L, the twists. A string that winds once around 
a twist cannot close, as it changes from dark to light or 
vice versa. That is, e charges become m charges if they 
are transported around a twist. Because of these global 
properties, it is clear that twists have a topological na- 
ture, in the sense that isolated twists cannot be created or 
destroyed locally. In particular, the parity of the number 
of twists in a region cannot be changed without altering 
the geometry in the boundary of the region. 

With twists, the labeling of e and m charges becomes 
inconsistent. But we can still fix a local labeling at a 
basepoint and make it global via a fixed set of paths. 
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FIG. 3: Several string operators (solid lines) and twists 
(crosses), in a top view of the system. The ordering of seg- 
ments in the string dictates which is drawn over which. Dot- 
ted lines are fermion string operators, two parallel string oper- 
ators of different type, (a) A process that creates/annihilates 
an e. (b) Two inequivalent closed string operators winding 
around a given twist. We choose the one to the left to define 
a+ and cr_. (c) Three equivalent string operators. The first 
equation uses property (ii) of string operators, the second uses 
also (i). (d) The strings 7^,7^,7'',... that define Majorana 
operators. The stars mark their endpoints. (e) The effect on 
7-' strings of braiding adjacent twists as in Fig. [T]^a) . 



Also, from a topological perspective only the location 
of the twists matters. In particular, if plaquettes are 
shifted along a circle enclosing a region without twists 
the topological properties of the system are unchanged. 

The possibility of considering dislocated lattices for ^ 
was already pointed out by Kitaev in [§] , where he argued 
the appearance of topological degrees of freedom. We will 
see not only that this is true, but also that twists have 
interesting fusion and braiding rules. 

Generalized charges — Twists act as sources and 
sinks for e fermions. Indeed, an e fermion can split in a 
pair of e and m charges, and if one of them winds around 
a twist they can annihilate, see Fig. [3Ja). This suggests 
that we cannot distinguish four topological charges in 
a region with a twist, as 1 and e, and also e and m, 
become identified. String operators make this more pre- 
cise. In a region with an odd number of twists we can- 
not define the two enclosing strings needed to distinguish 
the four charges. Instead, there is a single string 7 that 
winds twice around the region, see Fig. EJb). Since 7 
self-crosses, the ordering of its segments is relevant and 
we choose one of two topologically distinct orientations, 
see the figure. The eigenvalues of S*^ distinguish two new 
topological charges: cr+ for S-y = i and cr_ for S^ — —i. 
The pentagonal plaquette operator of Fig. HJb) intro- 
duces an energy difference between them. 

We have six generalized topological charges in the 
model: 1, e, m, e, cr+ and a-. These are not the charges 
of an anyon model, because anyonic charges are not mod- 
ified through braiding. Yet, they do have well defined 
fusion rules that can be recovered using string operators: 

cr± X (T± = 1 + e, (T± X fj^ = e + TO, 
CT± X e =CT±, (T± X e = cr-t X TO = (7^. (6) 

To exemplify the computation of these rules, consider the 
fusion (T_|_ X Let 7^, 7^ be the strings that wind twice 
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around each of the twists, as in Fig.[3jc), and let 7^, 7^ 
be two nonequivalent strings that enclose the two twists. 
Then, as demonstrated in the figure, in the absence of 
other excitations S^oS^l — —S^iSy2 = 1, which implies 
that the total charge can only be 1 or e. Since we can 
switch the charge by moving an external e "into" any of 
the two twists, the rule follows. 

Ising anyons — In view of (|6]), it is tempting to con- 
sider the subset of charges {l,(T+,e}: it is closed under 
fusion — with rules ([T]) via (T+ — > cr, e — — and a 
charge e remains the same after being transported around 
cr+. Could this subset of charges be regarded as an 
anyon model? Not really: we cannot attach an invari- 
ant meaning to Rea+Ra-^e, because we cannot directly 
compare 27| the path of an e that surrounds a (T+ and 
one that surrounds a trivial charge: the geometry is dif- 
ferent. Nevertheless, as we will show now, the braiding 
and fusion of cr+ twists reproduces that of a Ising anyons. 

If braiding twists is to make sense, we need to move 
them. For example, by adiabatically transforming the 
"geometry" of the Hamiltonian. In a different context, 
that of toric codes 0, [2lj , there is no Hamiltonian but 
just an encoding subspace that corresponds to the ground 
subspace of ([5]). There we can also change the geometry 
through code deformations 21, 2^ 2^. Although the ge- 
ometry of the Hamiltonian or code depends necessarily to 
a certain extent on the braiding history, at a topological 
level we can ignore this as long as different processes are 
not compared. Then, whatever the setting, with or with- 
out a Hamiltonian, we can analyze the effect of braiding 
following the evolution of closed string operators, which 
are topologically protected. We will focus on the braid- 
ing of (7+ twists, but the treatment of (7+ and cr_ twists 
together is very similar. 

The key to understand braiding is to introduce strings 
7^ such that the operators Cj := S^j realize the Majorana 
operators in ([2]). A suitable choice for the 'y^ is given in 
Fig. |3l[d). All these strings have the same endpoints, so 
that any product of an even number of Cj-s corresponds 
to a closed string configuration and gives a constant of 
motion. This is true in particular for the operators Cj :— 
—iCjCj^i, which give the total charge of the j-th and 
{j + l)-th twists. Namely, Cj = 1 if the charge is 1 



and Ci 



-1 if it is e, exactly as desired. Moreover, 



under braiding the evolution of the constants of motion is 
dictated by ([2]), as can be computed graphically starting 
from Fig. ([3|)(e). In conclusion, (t+ twists exactly mimic 
Ising anyons under braiding and fusion. That is, as long 
as process dependent phases are not involved, which fails 
for example in an interference experiment. 

Discussion — Beyond the particular example ad- 
dressed here, it would be interesting to explore twists in 
general anyon models. E.g., can twist braiding and fusion 
be computationally universal in a non-universal anyon 
model? Also, we have adopted an approach where twists 
are directly engineered, but they could appear randomly 



in systems where topological order emerges naturally. In 
the context of topological codes, the ideas presented here 
are tools to design and manipulate codes. E.g., they al- 
low to get planar versions of the topological subsystem 
codes introduced in (24} . 
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